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Classical chemical kinetics use rate-equation models to describe how a reaction proceeds in time. 

Such models are sufficient for describing state transitions in a reaction where coherences between 
different states do not arise, or in other words, a reaction which contain only incoherent transitions. 

A prominent example reaction containing coherent transitions is the radical-pair model. The kinetics 
of such reactions is defined by the so-called reaction operator which determines the radical-pair state 
as a function of intermediate transition rates. We argue that the well-known concept of quantum 
walks from quantum information theory is a natural and apt framework for describing multisite 
chemical reactions. By composing Kraus maps that act only on two sites at a time, we show how 
the quantum-walk formalism can be applied to derive a reaction operator for the standard avian 
radical-pair reaction. Our reaction operator predicts a recombination dephasing rate consistent with 
recent experiments [J. Chem. Phys. 139, 234309 (2013)], in contrast to previous work by Jones and 
Hore [Chem. Phys. Lett. 488, 90 (2010)]. The standard radical-pair reaction has conventionally 
been described by either a normalised density operator incorporating both the radical pair and 
reaction products, or by a trace-decreasing density operator that considers only the radical pair. 

We demonstrate a density operator that is both normalised and refers only to radical-pair states. 
Generalisations to include additional dephasing processes and an arbitrary number of sites are also 
discussed. 

PACS numbers: 03.65.Yz, 03.67.-a, 82.30.-b, 05.40.Fb 


I. INTRODUCTION 

It is known that many animal and insect species are ca¬ 
pable of sensing extremely weak magnetic fields. Of par¬ 
ticular interest amongst biologists, chemists, and physi¬ 
cists is the problem of how migrating species of birds 
use the Earth’s magnetic field for navigation pQ. The 
mechanism granting birds the ability to use the geomag¬ 
netic field for guidance is known as the avian compass 
and there is now substantial evidence that certain species 
(e.g. the European robin) do indeed possess this com¬ 
pass. To date, the most promising model of the avian 
compass is known as the radical-pair mechanism EHS], 
a chemical reaction which takes place inside a photore¬ 
ceptor molecule in the bird’s eye known as cryptochrome 
HHiij- The radical-pair model has been a platform on 
which many interesting theoretical investigations have 
been carried out since it was first proposed as a can¬ 
didate explanation for the avian compass fl3l [14] . One 
interesting problem, which has also been the subject of 
recent debate, is the form of its reaction operator. We 
will use quantum walks, which is essentially an elaborate 
theory of Kraus maps pans], to shed some light on this 
topic. This illustrates that quantum walks is a suitable 
framework for describing coherent chemical kinetics. 

The description of radical-pair reactions has con¬ 


ventionally been that of Haberkorn’s from 1976 El- 
This approach is phenomenological and based on argu¬ 
ing which of two existing differential equations for the 
radical-pair state should be preferred. The first is pro¬ 
posed by Johnson and Merrifield [18], and Evans et al. 
PS- The second is due to Pedersen and Freed [20] . We 
shall refer to Refs. QUO!] collectively as the Johnson- 
Evans equation for mere convenience. Haberkorn chose 
the Johnson-Evans equation instead of the one due 
to Pedersen and Freed by showing that Pedersen and 
Freed’s equation leads to negative eigenvalues for the 
radical-pair state whereas the Johnson and Evans ver¬ 
sion does not. Both differential equations for the radial- 
pair state can now be seen to contain terms making up 
what is known as the Lindblad form of master equations 
[21] [22], though neither is actually in Lindblad form. 
Haberkorn’s solution was to simply consider what the 
system state would look like if propagated using the two 
proposed state derivatives. The nontriviality in distin¬ 
guishing the two types of state evolution lies in the fact 
that it is not immediately obvious how to interpret the 
proposed state derivatives as opposed to the state itself. 
In more general and formal language, Haberkorn made 
his argument by referring to the map JC(t) which takes a 
system state from p(0) to p{t ), rather than the map’s gen¬ 
erator £, which is related to the map by JC(t) = exp(£t). 
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Hence the generator hence defines the derivative of the 
state, i.e. p(t) = Cp(t). In the context of chemical ki¬ 
netics, the superoperator C is referred to as the reaction 
operator. Of course, today, the Lindblad form is well un¬ 
derstood (at least in the quantum optics and quantum 
information community), so the problem pointed out by 
Haberkorn may not be perceived by some to be signifi¬ 
cant. However, there is still something to be learnt from 
Haberkorn’s work, which is that when considering non¬ 
unitary evolution the map JC(t) may be a more intuitive 
quantity to consider than the state derivative C because 
p(t) is given explicitly by JC(t) but only implicitly by C. 
One could in fact argue that this is also the reason why 
Lindblad’s result on the form of master equations is non¬ 
trivial EH [22]. This motivates us to use /C(t) instead of 
C in this paper. 

Haberkorn’s preferred reaction operator went unchal¬ 
lenged until recently [22], bringing the debate about 
its form into the limelight again [231130] . This has re¬ 
sulted in one side arguing in defence of Haberkorn and is 
now referred to as the conventional, phenomenological, 
or Haberkorn approach to radical-pair reactions j24H26] . 
A separate camp, called the quantum-measurement ap¬ 
proach to radical pairs has proposed two new reaction 
operators, due separately to Kominis [231 and Jones and 
Hore m- Of particular interest to us is the paper by 
Jones and Hore [27 who derived their reaction operator 
using Kraus maps m- The Jones-Hore equation pre¬ 
dicts a different singlet-triplet dephasing rate to the con¬ 
ventional approach of Haberkorn’s and has been the sub¬ 
ject of a recent experiment aimed at discriminating the 
two models |3T. This experiment showed the Jones-Hore 
equation to be inconsistent with the measured dephasing 
rate. In this paper we also use Kraus maps to describe 
the radical-pair kinetics but we obtain a dephasing rate 
that is consistent with Ref. m- a key factor in our ap¬ 
proach is the recognition that any intermediate transition 
in a multistate reaction involve only two states at once. 
Although this seems trivial, it is what separates our pa¬ 
per from the work of Jones and Hore because it implies 
that one can derive the map for a multistate reaction by 
composing two-state maps only. Maps for multistate re¬ 
actions derived in this way will thus be correct by virtue 
of the method (provided that we have the correct two- 
state maps). Two-state transitions are particularly well 
understood in quantum information theory since qubits, 
which are essentially two-state systems, are the central 
object of study. The toolbox provided by quantum infor¬ 
mation theory allows us to construct maps for multistate 
reactions that are robust to modelling errors. 

Our approach to the radical-pair reaction kinetics 
views the reaction as simply a system evolving between 
a discrete set of states in a probabilistic manner. Since 
such systems are analogous to random walks, our ap¬ 
proach to the problem is one of quantum walks. We re¬ 


view quantum walks below and point out the sense in 
which our version of quantum walks differ from those in 
the quantum-walk literature. 

Quantum walk was first introduced in the early 90s by 
Aharanov and coworkers [32] . They sought to generalise 
the idea of a classical walker who can only move left or 
right in discrete units along one spatial dimension to the 
quantum case. They managed to define a quantum walk 
as the analogue of a classical random walk by correlat¬ 
ing the system’s spatial coordinate to its internal degree 
of freedom such as spin, generically called a coin. The 
coin’s ability to be in a superposition of states can be 
seen to give rise to quantum walks although the use of a 
coin is actually not necessary [33] . Since then quantum 
walks have proven to be useful in quantum information 
[34] where they have found a variety of algorithmic ap¬ 
plications [35] in hitting [36] [37] and searching [38] [39] . 
Quantum walks were first introduced for closed systems 
which follow unitary evolution, but have recently been 
generalised to open systems which follow nonunitary evo¬ 
lution, called open quantum walks miii]. Such evolu¬ 
tion is described by a map that is completely positive 
and trace preserving m, and like their unitary counter¬ 
part, the maps defined in Refs. mm include changes 
in the internal degrees of freedom of the open system. To 
model the radical-pair reaction we propose an evolution 
map which makes no reference to any internal degrees 
of freedom. In this sense our model of the radical-pair 
reaction is similar to Ref. |33j with the exception that we 
allow for nonunit ary evolution. 

The rest of the paper is organised as follows. In Sec.[H 
we define the standard radical-pair reaction and review 
the different reaction operators that have so far been pro¬ 
posed. These results form the backdrop against which 
our approach to radical-pair kinetics should be consid¬ 
ered. We then introduce Kraus maps in Sec. Ill which 
writes maps in a particular form known as the operator- 
sum representation m Here maps describing processes 
known as amplitude-damping, dephasing, and unitary 
evolution are explained. We give a detailed exposition of 
the amplitude-damping map in Appendix [A] to illustrate 
how the operator-sum representation provides insight to 
the process which would not have come by so easily if 
the process was described using a Lindblad-form master 
equation. This forms our toolbox for describing radical- 
pair kinetics and is used in Secs. [IV] and [V] to derive a 
reaction operator for the standard radical-pair reaction. 
In Sec. [IV] we focus on reaction operators which can be 
described with only amplitude-damping maps. This cor¬ 
responds to the reaction operators reviewed in Sec. [II 
where only recombination processes are assumed to oc¬ 
cur. An interesting result here is the radical-pair density 
operator obtained from a partial trace over the chemical 
products. This gives a new state in the radical-pair basis 
which has been overlooked in previous models. In Sec.|V| 
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sufficiently far away in energy so that they can be safely 
neglected. We label the singlet and triplet product states 
as |P S ) and |P T ). The |S) —► |P S ) and |T) —► |P T ) 
transitions are then characterised by the respective rates 
ks and for- We have circled the singlet and triplet states 
in Fig. [I] with a dashed line to emphasise that the sys¬ 
tem comprises of only the |S) and |T) states. All known 
reaction operators are expressed in terms of these states 
so it will be convenient to define 


Fig. 1: The standard radical-pair reaction. The spin state of the 
radical pair oscillates coherently between singlet and triplet states 
under the hyperfine interaction with neighbouring nuclei spins. 
This oscillation can be modulated by an external magnetic field 
and shown to be sensitive to the magnetic field’s direction. This 
effectively modulates the amount of time the radical pair spends 
in the singlet state versus the triplet. Since each spin state decays 
to their own unique product, information about the magnetic field 
direction is then encoded in the concentrations of the singlet and 
triplet products. 


we generalise the results of Sec. [IV] to include dephas¬ 
ing and unitary dynamics and comment on the relation 
between our quantum-walk approach and previous work. 
Finally we summarise our results in Sec. |VI| and discuss 
its connection with a sequel paper and other relevant lit¬ 
erature. 


II. REACTION OPERATORS IN THE 
LITERATURE 

A. The standard radical-pair reaction 


Qs = |S)(S| , Qt = |T)(T| . (1) 

The green dashed boundary in Fig. [l]then serves to re¬ 
mind us that in the present discussion, Qs and Qt sums 
to the identity, i.e. Qs + Qt = 1. This may be useful to 
keep in mind for later as our approach extends the sys¬ 
tem Hilbert space to include the product states so that 
|S) and |T) no longer form a complete set. 

The problem is to determine the appropriate reaction 
operator which describes the changes brought upon p(t) 
due solely to the spin-selective recombination taking the 
spin states to product states. In particular, this recom¬ 
bination contributes to the singlet-triplet dephasing and 
we would like to determine what exactly this contribu¬ 
tion is. All other effects such as spin relaxation or ef¬ 
fects of molecular motion are ignored. Spin-dependent 
interactions such as the Zeeman, hyperfine, exchange, or 
dipolar are also ignored. Including additional processes 
other than the spin-selective recombination will tend to 
increase the dephasing rate. This will be assumed in all 
of the reaction operators which are reviewed next, which 
all describe Fig. [l] but without the coherent evolution. 


The body of literature discussed here refers to the stan¬ 
dard radical-pair reaction shown in Fig. [l] This is often 
referred to as a spin-selective recombination reaction |6]. 
It is spin selective because the reaction product depends 
on the spin state of the reactants, i.e. the radical pair, 
while the term recombination refers to the physical pro¬ 
cess by which the chemical products are obtained: The 
radical pair is typically formed through the transfer of an 
electron from one molecule, called the donor, to another 
molecule, the acceptor, creating a spatially separated pair 
of entangled spins. The formation of the chemical prod¬ 
ucts usually involve a back-transfer of this electron from 
the acceptor to the donor—a recombination of the elec¬ 
tron with the vacancy left on the donor molecule (see 
e.g. Refs. [6j 42 , 43 j for a more complete account of the 
radical-pair mechanism). 

In line with previous works pH [23 EE] we assume a 
minimal basis {|S), |T)} for the reaction where |S) de¬ 
notes the singlet state and |T) the triplet state with zero 
total spin. Physically this corresponds to the high-field 
limit where the triplet states with nonzero total spin are 


B. Proposed reaction operators of the standard 
radical-pair reaction 

1. Theoretical models 


Beginning with Haberkorn, the reaction operator pre¬ 
ferred by his positivity-preserving argument for the 
radical-pair state is given by 


pit) = £ h pit) 

~\ k T 


Qs pit) + pit) Qs 
Qt pit) + p{t) Qt 


( 2 ) 


Note that this is a non-trace-preserving equation for p(t), 
the effect of which is to describe leakage of the singlet 
and triplet populations over time at the rates ks and &t 
respectively. Equation Q also predicts a damping of the 
singlet-triplet coherence at the rate of (ks + for)/2. This 
can be seen by calculating the time derivative of the off- 
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diagonal elements of p(t): 

Pst(£) = — ^(^s + ^t)pst(^) , (3) 

where Pst(^) = (S|p(t)|T). This reaction operator was 
first questioned by Kominis [23] who argued that the 
radical-pair reaction is analogous to two coupled quan¬ 
tum dots under continuous measurement by a point con¬ 
tact. Kominis then derived a reaction operator using 
similar methods as Refs. [44H46] which gives 

p(t) = C K p(t) 

= £h p(t) + ks Qs pit) Qs + kr Qt p(t) Qt ■ (4) 

We have written Ck in terms of Ck to emphasise that the 
difference between 0 and © lies in the terms Qs p(t) Qs 
and Qt p(t) Qt- The addition of these terms puts Komi- 
nis’s result in Linblad form making the evolution trace 
preserving. This means that Ck does not describe the 
loss of singlet or triplet populations as in (J2|. Kominis 
therefore augments his description of the radical-pair ki¬ 
netics by an additional equation in which the radical-pair 
population is given by 

N{t + dt) = N(t) [l - ps(t) - pr(t)] , (5) 


where 


ps(t) = k s Tv[Qs p(t)]dt , PT = k T Tr[Q T p(t)]dt (6) 


are the respective probabilities of a transition from either 
|S) to |P S ), or |T) to |Pt) in the infinitesimal interval dt. 
Equation Q does however predict the same dephasing 
rate as (J2|. This is obvious from © since Qs p(t) Qs and 
Qt p(t) Qt do not contribute anything to (S|p(t)|T) due 
to the orthogonality of |S) and |T). 

Spurred on by the measurement analogy made by Ko¬ 
minis, Jones and Hore m attempted a derivation of 
the radical-pair reaction operator using the operator-sum 


representation (see Sec. III). Their result is given by 


to be too small a difference to be detectable in an ex¬ 
periment EZ3, Maeda and colleagues have recently man¬ 
aged to place the two models under experimental scrutiny 
m using pulsed electron paramagnetic resonance spec¬ 
troscopy m- We briefly review some key elements of 
this experiment below. 


2. Experimental discrimination 

The experiment reported in Ref. m uses radical 
pairs in a modified version of the carotenoid-porphyrin- 
fullerene triad molecule of Ref. [4]. This system has 
previously been demonstrated to exhibit the anisotropic 
chemical response to Earth-strength magnetic fields re¬ 
quired for the avian compass H HHJ. This model sys¬ 
tem also minimises the singlet-triplet dephasing due to 
processes other than recombination and has a &t much 
smaller than ks for some temperatures (between 200- 
240 K). This means the dephasing rate in this tempera¬ 
ture regime is approximately ks/2 according to £h, and 
ks according to £jh- The rate ks was then measured 
for this temperature range. It was shown that for cer¬ 
tain temperatures (a range of about 30 K), the dephasing 
rate from the Jones-Hore model lied above the measured 
values while the rate from Haberkorn’s model always re¬ 
mained below. In practice there will be other uncon¬ 
trollable processes which tend to increase the dephasing 
rate so the measured values of dephasing will not come 
solely from the |S) —> |Ps) recombination. This means 
that any reaction operator must produce a dephasing rate 
which lies below the measured values of ks for all temper¬ 
atures and therefore suggests the recombination kinetics 
according to £jh to be incorrect. In the next section we 
show how the idea of quantum walks can be used to de¬ 
rive a reaction operator with a dephasing rate consistent 
with the experimental data of Ref. [31 . 


III. OPERATOR-SUM REPRESENTATION 


p(t) = C m p(t) 

= — (ks + &T )p(t) 

+ k s Qt p(t) Qt + Qs p{t) Qs • (7) 

Similar to the phenomenological approach given by ©, 
this equation also does not preserve the trace of p(t). As 
already mentioned, this is attributed to the loss of the 
singlet and triplet populations to the reaction products. 
It can be seen that 0 gives this loss rate at ks and &t 
for the singlet and triplet states respectively as expected. 
However, the Jones-Hore reaction operator Cjk gives a 
dephasing rate which is the sum of the recombination 
rates, i.e. ks + &t> rather than the average predicted by 
£h (o r Ck)- Despite being perceived by Jones and Hore 


Kraus published his theory of general state changes in 
quantum mechanics m in which he asked what form 
must a superoperator Q take if it is to map a physically 
valid state p to another physically valid state p '. Note 
that we have not actually mentioned time so p' can be 
the state of a quantum system after some abstract op¬ 
eration, not necessarily a state at a particular instant in 
the future of p (although we will use it to propagate the 
system in time). The answer to the question just posed 
is that Qp must be of the Kraus form, also known as the 
operator-sum representation of Qp. Kraus’s result there¬ 
fore has the power to describe a large class of state tran¬ 
sitions without referring to the underlying physics which 
makes the theory operational. This is what gives the 
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Kraus formalism its versatility. To describe time evolu¬ 
tion we simply associate p with the system state at some 
arbitrary time t, and p' at some later time t + At. The 
operator-sum representation can then be stated as 

p{t + At) = S(A t)p(t) 

N 

= ££<"> (At) p(t)[£<"> (At)] f , ( 8 ) 

n= 1 

where the set of Kraus operators {K^ n \At)} n satisfy the 
condition 


N 

]T [. k {n) (At)] f k {n) (At) = i. (9) 

n= 1 

The theory also gives a prescription for calculating the 
probability that event n is observed, given by 

Pn (At) = Tr| [K^ (At)] f (At) p(t )| . (10) 

The sum in <© can be understood to be over states con¬ 
ditioned on events (indexed by n) that may be observed 
in an interval At, and hence its connection to measure¬ 
ments. Condition © is then equivalent to the conserva¬ 
tion of probability expressed by 


N 


)T Pn(At) = 1 


( 11 ) 


n= 1 


Note that (10) is simply the norm of the post¬ 
measurement state so that (J 8 | may also be written as 
an average over normalised conditioned states: 


N 

p{t + A t) = J2 P (n) ( A 4 P (n) (t + A *) > (12) 

n= 1 


where 

p (n) (t + At) = K (n) (At)p(t ) [ K (n) (At)\ f /p (n) (At ) . 

(13) 


A. Amplitude damping 

The first idea that we will borrow from quantum infor¬ 
mation theory to describe the radical-pair reaction is the 
amplitude-damping Kraus map. A derivation of this map 
can be found in Ref. m but is given in terms of a pho¬ 
ton incident on a beam splitter. This is an example of a 
process that can be described by the amplitude-damping 
map but we believe an explanation involving only ideas 
from probability theory and basic quantum mechanics 
to be more direct and fitting for developing the reaction 
operator. We have thus included such a detailed expo¬ 
sition in Appendix [A] and provide only a sketch of the 
amplitude-damping map below. 


721 

• -► • 

IV’i) 1^2) 


M 21 

A 

• • 

IV’i) 1^2) 


(a) 


(b) 


21 


IV’i) 


1 ^ 2 } 


(c) 


Fig. 2: (a) Amplitude damping from |^i) to \ip2)- This map trans¬ 
fers population from one state to another, (b) Dephasing between 
|'0i> and |b2>. The map tends to localize the random walker onto 
l^i) or | / 02) (i.e. turn a superposition of I'lpi) and ^ 2 ) into a mix¬ 
ture). The triangle can be thought of as a wedge driven into a 
line that connects the two states being dephased. (c) Coherent os¬ 
cillations between states |'0i) and \ip2)‘ The interconversion rate 
between states |'0i) and ^ 2 ) is given by 2 £21 [see ( |42| )]. 


Assuming first for simplicity that we have a two- 
dimensional system described by |^i) and IV^)- Aside 
from its Hilbert space dimension, the system is other¬ 
wise general, and the states | 7 / 7 ) and ^ 2 ) are arbitrary. 
We now wish to describe the change in the system state 
over some time interval At, say from t to t-\- At, allowing 
for the possibility that a transition from to ^ 2 ) ma y 
occur during this interval. This is schematically shown 
in Fig. [2] (a). 

We can think of Fig. [ 2 ] (a) as describing a possible 
change of state for a single molecule out of an ensem¬ 
ble of n identical molecules. Thus some fraction of the 
molecules will jump from state |Vq) to ^ 2 ) over some 
finite time interval. By regarding the transition as a ran¬ 
dom event we can characterise the process by a prob¬ 
ability, 721 (At), for a given molecule to go from state 
| vpi) to 1 ^ 2 ) in time At. If we prepare all n molecules in 
state |Vq) and observe them for a time of At, then the 
fraction of molecules that make the transition to state 
|^ 2 ) is given by 721 (At). We should then expect that 
the longer we observe the process the greater the num¬ 
ber of molecules to jump to state IV^)- In the long-time 
limit, all n molecules end up in state 1 ^ 2 ) so we expect 
721 (At) —)► 1 as At —> 00 . Conversely if the process is 
only observed for a very short interval then we would not 
expect many molecules to have jumped to state 1 ^ 2 )• We 
thus expect 721 (At) —> 0 for At —> 0. We denote the 
map describing such a process by M. 2 i{At). It has the 
operator-sum representation with two Kraus operators: 

p{t + At) = A4 2 i(At) p(t) 

= M 2 ( 1 1 ) (Af)p(f)[M 2 ( 1 1 ) (Af)] t 







6 


+ Mg\At) p(t)[M£\At)]\ (14) 

where 

^21 ( At ) = \/72i(A t) |^2) (^i | , (15) 

M^(At) = |- 02 > <- 02 | + \A ~ 721 (A t) \i>i)(ipi\ • (16) 


Its effect on an arbitrary state can be seen most directly 
by calculating the matrix representation of (14) in the 
basis { 1 ^ 1 ) 5 1^2) }• This gives the 2x2 matrix 


p(t + A t) 

= (Pn(t) -72i(Ai)/j n (i) y/l - 721 (At) p 12 (t) \ 

\ \/l — 721 (■ At) P 21 (t) p 22 (t) + 721 (At) Pll(t) J 

(17) 


where we have defined pjk = (^j\p\ipk)- The population 
transfer from state |Vq) to ^ 2 ) is apparent on the di¬ 
agonal terms [54 in A fraction 721 (At) has been 

subtracted from pn(t) and added to P22(t). Note that 
the off-diagonal terms of p(t) are also affected by this 
process. Unless 721 (At) = 0 the population transfer will 
simultaneously reduce the coherence between |?/q) and 
|^ 2 )- This can be seen from AppendixjA] where we argued 
about the form of (15) and (fl 6 |) without ever referring 


to the system coherences. The decay of the off-diagonal 
elements in © should thus be taken as a consequence 
of the population transfer. This is a crucial difference be¬ 
tween our formulation of the reaction operator and that 
of Ref. [27] where the decay of coherences was put into 
the system evolution by hand. When l^i) represents a 
state of higher energy than IV^), © is said to describe 
a dissipative process (hence the name amplitude damp¬ 
ing) . In this case © captures the well-known idea from 
open-systems theory that dissipation implies decoherence 


We now generalise the amplitude-damping map to a 
system with N states. By essentially the same argument 
as in Appendix [A] the map describing a transition from 
state I'lpk) to | x/jj) is simply given by 

M jk (At) p{t) = M$iAt)p(t) [M${At)] ] 

+ M}l\At)p(t)[M$(Atj\\ (18) 

with the Kraus operators 

{At) = yj7 j/c(At) Qjk , (19) 

(A t) = P k + ^/l- 7 ift(A t) Q k , ( 20 ) 

where jjk(At) G [0,1] for all j, k , and At. For ease of 
writing we have defined 


Since {\^k)}k=i * s assumed to form a complete set for an 
arbitrary AT, the identity operator may be written as 


N 

ijv = Qk • 

k=1 


(23) 


Note the order of subscripts in (18)— ([20]) is important. 
Reversing the order of the subscripts reverses the direc¬ 
tion of the process. It can be verified that © is repro¬ 
duced by taking j = 2, k = 1, and N = 2 in (fl8])-(23). 


Although we have parameterised the amplitude-damping 
map Adi j (At) by the probability 7 ^ (At), it is usually the 
rate kij at which the transition occurs that is measured 


or estimated. It is also a more useful quantity to use 
when expressing the amplitude-damping map in differ¬ 
ential form. If we let kij be the fraction of particles that 
jump from | i/jj) to | ipi) in one second, then the fraction 
of particles that make the transition in time At is simply 
given by 


lij(At) = k^ At . 


(24) 


It will be useful to express (18) in differential form. 
The two-dimensional case studied in Appendix [A] can be 
generalised to the case of N states [see ( A15| )]. The dif¬ 
ferential form of Mij(t) is given by 


dp__ r 
dt " 


— kij 


P(t) 

Qij p(t) Qjj 


\ Qj p(P) ~ \ p(t) Qj 


(25) 


This can be written in the Lindblad form if preferred by 
noting that Qj = QjiQij = Q\jQij- We can formally 
express the amplitude-damping map as 

M ij {t) = e c ^ t , (26) 


where Cij is said to be the generator of Adij(t). The 
reader may proceed directly to Sec. [IV] at this point if 
he/she wishes as the results there only require knowl¬ 
edge of the amplitude-damping map. We will introduce 
the maps for dephasing and coherent evolution next but 
they will not appear until Sec. |V|and hence can be read 
then. Lastly, note that we have suppressed the depen¬ 
dence of superoperators (either maps or generators) on 
the transition rates (or probabilities). This is because 
here we are thinking of the transition rate as a fixed num¬ 
ber, a parameter which defines the system as opposed to 
time, which is an independent variable. In order to keep 
our notation simple we will continue to suppress param¬ 
eter dependencies unless otherwise stated. 


B. Dephasing 


Qjk \^j)(^k\ 5 
Qk = \^k)(^k \, Pk = i — Qk • 


( 21 ) 

( 22 ) 


We saw in © that population decay in one of the 
states led to decoherence. However, decoherence may 
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also occur without population decay and this is known 
as dephasing, or phase damping in analogy to amplitude 
damping (so called because it results from a loss of infor¬ 
mation about the relative phases between the different 
basis states | ip k ))- This process is represented by the 
symbol shown in Fig. |2](b). The system evolution under 
dephasing over a time of At can be represented simply 
by a 2 x 2 matrix for a two-state system as 

p{t + At) 

= V 2 i(A t) p(t) 

pn(t) y/l ~ P2i(At) p 12 (t) 

\/l - M2l(A£) P21 (t) P22 it) 

(27) 

where /x 2 i(A t) G [0,1]. Note that only the coherences 
(the off-diagonal terms in p) are damped. The generali¬ 
sation to a system with N states can be stated simply in 
Kraus form as 

V*fc(A t)p(t) = v£\a t)p(t) [V$\A i)] f 

+ V%\At)p(t)[V%\At)]\ (28) 

where 

V$\At) = yJn jk (At) Qk , (29) 

V$( At ) * h + y/l-p jk {At) Q k . (30) 

As with amplitude damping, we can work with the rate of 
dephasing rather than with pjkiAt). Denoting the rate 
of dephasing between states | !pj) and | ipk) as qj k , we can 
write 




PjkiAt) = q jk At . (31) 


The evolution under the dephasing map Vjk{t) can then 
be expressed by the differential equation 


It ~ Sjk p ^ 


— Qjk 


Qk pit) Q\ - - Qkpit) - - pit) Q k 


which defines the generator for Vjfc(t): 

Vjk(t) = e s > kt . 


(32) 


(33) 


C. Coherent evolution 


coherent evolution between two states graphically by us¬ 
ing a green two-way arrow as shown in Fig. [2] (c). In gen¬ 
eral, coherent oscillations between states \ifjj) and \^k) 
can be generated by a Hamiltonian of the form (for j ^ k) 

Hjk — Qj T Qk T ftjfc ( Qjk + Qkj ) 5 ( 34 ) 

where uj k = ('ipk\Hjk\'ipk) is the expectation value of Hj k 
in the state IVtfe)- The coupling between states l^j) and 
1 2pk) is denoted by Qj k . Note that for Hj k to be Hermi- 
tian Qjk must be real and symmetric with respect to its 
indices, i.e. 


Q jk = Q* k = ^kj • (35) 

Unitary evolution can be understood as a special case of 
the Kraus decomposition © with only one Kraus oper¬ 
ator: 


Ujk(At) = e-^ At , (36) 

where we have set h m 1 for convenience. The actual 
evolution over a time of At is then effected by the map 

U jk {At) pit) = U jk (At) p{t) u] k (At) . (37) 

The differential form of lAj k (t) has the familiar commu¬ 
tator form: 


= 7 ljkp(t) = -i\H jk ,p(t )] . 


dp 
dt 

We can also express ( |37| ) in terms of 1Zj k as 

lAjkit) 


'TZnkt 


(38) 


(39) 


Just as we parameterised the amplitude-damping and 
dephasing maps by their probability of occurrence, we 
can similarly parameterise unitary evolution by 


a jk (At) = \(ipj\Uj k (At)\if ) k )\ 2 . 


(40) 


This is the probability of making a transition to | t/jj) 
after time At assuming the system was initially in |^). 
If we wish to relate ajkiAt) to the transition rate under 
unitary evolution then an explicit expression for (40) is 
required. This can be shown to be 


ft? 


ajkiAt) — 2 [l cos i%Cjk At)] , 

^jk 


(41) 


where we have defined 


The map (28) describes pure decoherence and is use¬ 
ful for modelling processes which counter act any coher¬ 
ent evolution of the system that may occur in the ba¬ 
sis {\^k)}k=v The singlet-triplet inter conversion in the 
radical-pair mechanism is one such process. We depict 


(jk — ^ \Ji^k Uj ) 2 + 4 Q 2 k • (42) 

Equation © tells us that 2Q k can be defined as the fre¬ 
quency at which the system oscillates between states | ipj) 
and | ip k )- Note that this depends on both the coupling 












IV’i) 1^3) 

• • 


721 


743 


• • 

1^2} |V>4> 


Fig. 3: The standard radical-pair reaction without coherent evo¬ 
lution represented as a graph. Black arrows represent population 
transfer, which here is associated with the recombination reaction 
of the radical pair. 


between | ipj) and 1 ^) (i.e. Cljk) as well as their sepa¬ 
ration (given by \wk — ojj\). We can also see from p0| ) 
that increasing \ujj — Wk\ lowers the peak of the transition 
probability between | i/jj) and |^). The proof of (42) [and 
hence (40)] will be presented in a sequel paper where it 
is actually used to simulate an example quantum walk. 
Here it is sufficient to see how ajk is related to the rate 
of the process and the effect of varying and Cljk- 


IV. RADICAL-PAIR RECOMBINATION 
REACTION AS A QUANTUM WALK 

A. Reaction operator 


The quantum-walk formalism visualises state transi¬ 
tions in a quantum system as a network of nodes (repre¬ 
senting states) connected by edges (representing transi¬ 
tions), called graphs. Such models have a wide applica¬ 
bility because the nodes can represent abstract degrees 
of freedom, such as a spatial coordinate, or in our case, 
the state of a molecule. We therefore begin our quantum- 
walk model of the reaction outlined in Fig. [I] by simply 
representing the different radical-pair and product states 
as nodes on a graph labelled according to: 


IS) = IV’i), 

|T) = 1^3) , 

(43) 

IPs) = 1^2) , 

ip t > = m , 

(44) 

while the rates are taken as 



(N 

III 

m 

&T = &43 • 

(45) 


The states in (43) and (44) are assumed to represent dis¬ 
tinct stages of the radical-pair reaction, therefore we take 
{\^k)}k=i be an orthornormal basis. This is shown in 
Fig- HI As explained in Sec. [TTJ here we concentrate only 


on the recombination processes which take the system 
from the singlet to singlet product (|Vq) — > IV^)), and 
from the triplet to the triplet product (^ 3 ) —> ^ 4 )). 
Extensions to include dephasing and coherent evolution 
will be covered in Sec. El 

Our goal is to derive a reaction operator C. This means 
that we should express the evolution of p(t) in differential 
form. As before, we can simply obtain such an equation 
by propagating p(t) over an infinitesimal interval dt, ex¬ 
cept now there is more than one process happening at a 
time. This can easily be dealt with by using a single map 
/CQw(dt) composed from a series of maps. Each map in 
^qw(dt) represents a particular process in the system. 
Attributing each recombination process to an amplitude- 
damping map, we describe the evolution of p(t) by 

p(t + dt ) = K,Qw{dt) p(t) = M^{dt) M2i(dt) p(t) , 

(46) 


where on using (26) we have 


/Cqw (dt) = e Ci 3 dt e C2ldt = e {Ci3+C2l)dt . 


(47) 


Note the second equality follows because dt is infinitesi¬ 
mal. This can be seen by expanding the exponentials and 
neglecting terms on the order of dt 2 . This also means that 
reversing the order of exp (£43 dt) and exp(£ 2 i dt) does 
not change the second equality of ©. Since £21 and £43 
are independent of time, ( |47| ) shows that £43 + £21 is the 
generator of /Cqw(^) f° r any finite t [otherwise it is only 
the generator of }Cqw(dt)]. This is simply equivalent to 


= £qw/>W = (£21 + C^)p{t) . 


We can now apply (25) to (48) to obtain 


£qw p(t) 
= &21 


k 


Q21 p(t) Q 21 ~ 
Q 43 p(t) Qls 


43 


\ Qi p(t) ~ \ p(t) Qi 
~ \ Qa pit) - t p(t) Q 3 


(48) 


(49) 


A few remarks can be made about our result directly from 
the form of (49) but it will be easier to refer to its matrix 
representation in the site basis. This is a set of rate 
equations which includes the coherences between sites 
as well. We follow the standard convention of writing 
matrices where p mn = ('ipm\p\'ip n ) denotes the element at 
the ra th row and 77 th column of p. The matrix form of 
(49) can then be easily verfied to be given by 
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dp 

dt 


—^21 pn 
~ \ &21 P21 

~ \ (&21 + k^s)psi 
V \ &21 P41 


- \ k 21 Pl2 
&21 Pll 
" | &43 P32 
0 


(&21 + &43)/?13 “ | &21 Pl4 \ 


-1 &43 P23 
-&43 P33 
— i &43 P 43 


0 

- 2 &43 P34 
&43 P33 


(50) 


7 


r 


It is obvious from the matrix form of £qw p that the 
spin-selective recombination reduces all coherences ex¬ 
cept for the coherences of the two product states, given 


by P 24 and P 42 to be zero [recall (43) and (44)]. Of spe¬ 
cial importance is the decay of the singlet-triplet coher¬ 


ence. This is given by pi 3 in (50) which can be seen to 


be consistent with the experiment of Ref. m- Note that 
because p is Hermitian it follows that p is also Hermitian, 
so referring to pi 3 is the same as referring to P 31 . 


Adding the diagonal elements of (50) we see that (49) 


is trace preserving. In particular we see that 


dpn dp22 


dt 


dt 


= &21 Pll , — 


dp33 ^P44 


dt 


dt 


= &43 p33 - 

(51) 


That is to say, the rate at which singlet-state radical 
pairs are lost due to recombination is exactly balanced 
by the rate of increase of the singlet product. Recall from 
Sec. |IIB that previous treatments on the radical-pair ki¬ 
netics use trace-decreasing reaction operators which refer 
only to l^i) and ^ 3 ) [see © and ©]• It has been noted 
in Ref. {2? that such reaction operators produce singlet 
and triplet populations which satisfy (51) and thus poses 
no problem. However, a description in the minimal ba¬ 
sis { |Vh) 5 1^3)} still fails to account for coherences be¬ 
tween the radical pair and product states. Jones and 
Hore have argued that such coherent superpositions be¬ 
tween the radical pair and products decohere very quickly 
and is therefore consistent with a model in which they 
are neglected m- However, on accepting (|5Q|), we see 


that coherences between the radical pair and products in 
fact decay at a rate less than the radical-pair dephasing 
(e.g. P 21 compared with P 31 ) so the remark by Jones and 
Hore is not actually correct. Nevertheless, a model in 
which the product states are neglected is still permissi¬ 
ble so long as the radical-pair populations and coherences 
do not depend on the populations and coherences of the 
products. This is clearly true from the matrix form of 
£qw so we can write down such a reaction operator di¬ 
rectly. For ease of comparison with previous results we 


express this reaction operator in the notation of Sec. IIB 
[recall ©, @ and @]: 


£qw pit) = pss{t ) |S)(S| + psrit) |S)(T| 

+ PTs(t) |T)(S| + |T)(T| . (52) 

We have used an overbar on £qw to indicate that it is 


no longer trace preserving. We can simply read off pss. 


Ptt, and pst from (50) to get 


£qw p(t) = - ksQs p(t) Qs ~ 2 {ks + &t) Qs p(t) Qt 

— x (ks + k T ) Qt pit) Qs — k T Qt p(t) Qt- 

(53) 

The reader may have already noticed that pss 5 Ptt, and 
Pst are in fact the same as those given by the Haberkorn 
reaction operator £h given in which means that ( [53] ) 
and Q should in fact be the same. This can be shown 
by collecting terms porportional to ks as one group and 


terms proportional to fcx as one group in (53): 


£qw p(t) = - ks 

1 


Qs p(t) Qs + 2 Qs P(t) Qt 


Qt p{t) Qs 


- k 3 


Qt pit) Qt 


2 Qt p{t) Qs + 2 Qs P(t) Qt] 
Qs pit) + p(t) Qsj 


~2 ks 
“ 2 kT 


Qt p(t) + p(t) Qt 


(54) 


(55) 


The second equality is obtained by using Qt = I 2 — Qs 
and Qs = I 2 — Qt in the terms proportional to ks and 


kT respectively in (54). Note the identity operator carries 
a subscript 2 because it is now only an identity on the 
subspace spanned by the singlet and triplet states. The 
resolution of the full identity operator, I 4 , requires all 
four states of the radical pair and products. This is why 
we have circled the singlet and triplet states in Fig. [l] in 
Sec. |II A| We have thus derived the conventional spin- 


selective recombination operator using the operational 
and systematic treatment of quantum walks. 


Writing (53) in the form of (54) also makes the compar¬ 
ison with the Jones-Hore reaction operator easier. This 
is because, as argued by Jones and Hore, an alternative 
derivation of their result given by 0 begins with 

£jh p(t) 

= -k s Qs p(t) Qs + Qs p(t) Qt + Qt p(t) Qs 
— kT Qt pit) Qt + Qt pit) Qs + Qs pit) Qt - (56) 
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The difference between (54) and (56) can thus be seen 
in the coefficient of the cross terms [i.e. Qs p(t) Qt and 
Q t p(t) Qs]- The singlet-triplet dephasing rate was thus 
incorrectly posited at the outset in their argument. Al¬ 
ternatively, this can also be attributed to an incorrect 
formulation of the Kraus operators in the minimal ba¬ 
sis. Our approach on the other hand begins with all four 
states of the radical pair and products. This allows us 
to focus on describing the |S) —> |Ps) and |T) — > |Pt) 
transitions with the dephasing between |S) and |T) oc¬ 
curring as a consequence. Since the |T) —> |Pt) and 
|S) —> |P S ) transitions are identical we in fact only have 
to find the correct operator-sum representation for one 
of them and apply it twice to p(t) to obtain the map for 
the full reaction as explained in (|46|)-(48). This makes 
our quantum-walk approach less prone to modelling er¬ 
rors. Next we illustrate this point further by using the 
quantum-walk idea to obtain the appropriate map in the 
minimal basis by ignoring the chemical products. 


B. Ignoring chemical products 


If we are interested in deriving the Kraus map for only 
the radical-pair state then we need some way of captur¬ 
ing the effect of the decay from the radical pair to the 
products but without including the products in p(t). The 
above treatment of first deriving (50) and then reading 


off the equations of motion for the radical pair provides 
one way of achieving this result. Here we show an al¬ 
ternative method which is also based on quantum walks. 
The approach here is to find the operator-sum representa¬ 
tion for the radical-pair state with the chemical products 
ignored. 

Although it may sound contradictory, we will begin by 
including the product states in p(t). However, we will 
lump the states |Pg) and |Pt) into a single state which 
we denote as |P). This is shown in Fig.|4| The singlet and 
triplet states are defined as before in (43), but instead of 


(44) and (45) we now have 


|p> = IV*), 

ks = &21 5 &T = &23 


(57) 

(58) 


The reason for introducing Fig. [4] is twofold: First, it is 
easier to consider Kraus operators for a trace-preserving 


map. The evolution corresponding to (55) can then be 


extracted by using just one of the Kraus operators and 
hence will be trace decreasing. For this purpose it is 
sufficient to introduce only one additional state. This 
will then allow us to use the quantum-walk approach by 
composing two trace-preserving maps corresponding to 
the transitions shown in Fig. [4j The second reason for 
introducing Fig. [4] has to do with what is meant by “ig¬ 
noring the products”. In the language of quantum prob¬ 
ability, any unobserved degrees of freedom in a system 


IV’i) 




721 



1 ^ 3 } 



723 


1 ^ 2 ) 


Fig. 4: A simplified graph sufficient for deriving the Kraus op¬ 
erator corresponding to a trace-decreasing reaction operator. Due 
to the decoupling of the matrix elements in ( |88| ), this graph can in 
fact be used to simulate the partial trace provided that we ignore 
the coherences between the product and radical pair and that we 
understand the distinction between |P) and \ni = 0, 77,3 = 0) [see 
Sec. |IV B 2| especially the discussion between ( |88| ) and ( |89| )]. 


can be traced over to give a so-called reduced state. This 
procedure is called a partial trace and is known to be a 
trace-preserving operation. We will find that this gives a 
reduced state for the radical pair that seems equivalent 
to Fig. [4] but is in fact subtly different. It will thus be 
convenient to refer to Fig. [4] for comparison. We discuss 
these two ideas below. 


1. Evolution from discarding products 


The operator-sum representation of (55) can be de¬ 


rived by considering the quantum walk shown in Fig. [4j 
which can easily be described by 

p(t + dt ) = M(dt) p(t) = M23{dt) M2i(dt) p(t) . (59) 

It can be shown in general that a composition of two 
Kraus maps is another Kraus map. This means that 
we can define the Kraus operators for the total map 
M(dt) by substituting the operator-sum representation 


of M. 23 {dt) and M. 21 (dt) into (59). This gives us 


M(dt) p(t) = Y M^ n \dt)p(t) [M^ n \dt)Y , (60) 

n=0 

where we have defined 

Ah°) (dt) = Af$ (dt) M 2 ( f (dt) , (61) 

M (1) (dt) = M 2 ( 3 } (dt) (dt) , (62) 

M (2) (dt) = M 2 ( f (dt) M 2 ( f (dt) , (63) 

M (3) (dt) = M 2 ( f (dt) (dt) . (64) 

Using © and ( |20| ) we find that 

M^(dt) = 0, (65) 

which means that it is redundant and we need only three 
Kraus operators to describe Fig. [4j Using (24) and the 
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binomial expansion to first order in dt in (62)—(64) [see 
(A12)] we find 


M^\dt) = \f k 2 \ dt Q21 , (66) 

MW(dt) = y/fadtQ 2 3, (67) 

(dt) = i - ^ Qi dt - Q 3 dt . (68) 


We should note that (66)—(68) could have also been ob¬ 


tained by generalising the argument in Appendix [A] to 
three states. In this case we would actually just write 
down three Kraus operators directly. It is simple to check 


that (66)—(68) satisfy the condition to be a valid set of 
Kraus operators [see (| 8 |)-(fT 3 |) from Sec. IIIIl: 


[M (n) (dt)] + M {n) (dt) = i . 


(69) 


n= 1 


The advantage of using (60) is that it decomposes p{t-\- 
dt) into a sum of states 


p (n) (t + dt) = M< n > (dt) p(t) [M (n) (dt )] f , (70) 


conditioned on the outcome n of a measurement per¬ 
formed at time t. The measurement is devised to give us 
information about the transitions occurring in the sys¬ 
tem so that we associate the n = 1 outcome with the 
\ipi) —> 1 ^ 2 ) transition, n = 2 with 1 ^ 3 ) —> |^ 2 )j and 
n = 3 with no transitions. It is easy to see that ( [66] ) 
describes a jump from state |Vq) to |^ 2 ) giving the con¬ 
ditional state 

p^ (■t + dt) = k 2 i Q 21 p(t) Q 21 dt 

= l2i{dt) pn(t) ^ 2 )(^ 2 1 • (71) 

Note that 721 (dt) pn(t) is the trace of p^\t + dt) so 
that upon normalisation the evolved state is simply |^ 2 )- 
Similarly applying ( |67| ) we get 

p (2) (t + dt) = k 23 Q 23 p(t) Q 23 dt 

= 723 (dt) p 33 (t) 1 ^ 2 ) ( 1 P 2 1 ! (72) 


which is seen to describe a jump from state I'lps) to |^j 2 )- 
This leaves (68) which gives 


p {3 \t + dt) = p(t) - t k 2 1 Qi p(t) + p(t) Qi 


dt 


~ 2 h3 


Q 3 p(t) + p(t) Q 3 dt. ( 73 ) 


Note that this is the evolution given by Haberkorn’s equa¬ 
tion. The effect of (73) can be seen directly by applying 
it to different initial states. Letting p(t) = |Vq)(Vh| we 
find 


(74) 


where 1 — k 2 i dt is the trace of p^ (t+dt). A similar result 
can also be seen by letting p(t) = \^ps)(^> 3 1- If the system 
is in the product state then we expect it to remain in the 
product state forever since there is no process to take 
the system out of \^ 2 ). Indeed, setting p(t) = |^ 2 ) (V^l 
we find 


p (3) (t + dt) = \tp 2 ) (V>21 • 


(75) 


The evolution described by ([73]) is thus conditioned on 
the absence of recombinations. That ( [68] ) describes 
radical-pair evolution conditioned on no recombinations 
is not surprising in view of Appendix [A] where it can 
be seen to be so by construction. The Kraus operator 
consistent with the conventional description of radical- 


pair kinetics in the minimal basis is thus given by (68) 


and is the measurement approach that Jones and Hore 
sought after in Ref. [27]. That such an evolution equation 
is given by conditioning on unrecombined radical pairs 
has also been noted in Ref. m, but with the operator- 
sum formalism incorrectly applied as the experiment of 
Ref. [31] has now shown. 


2. Evolution from tracing out products 


Here we would like to derive an equation of motion 
for the radical pair from tracing over the products in 


(49). The partial trace is a formal procedure for ob¬ 


taining a density operator with the products ignored 
and is quite different to simply discarding the prod¬ 
ucts. Note that the partial trace is not simply the sum 
(^ 2 \p(t)\^ 2 ) + (^ 4 |p(t)|^ 4 ), but rather an operation de¬ 
fined only on systems with a tensor product structure. 
As such, we require that the Hilbert space H of p(t) in 


(49) be of the form 


HI = Hr (£) Hp , (76) 

where Hr is the Hilbert space for the radical pair and 
Hp is the space for the products. We can then derive the 
state of the radical pair with the products ignored from 

p R (t) = Tr P [p(t)] = ^ <: Xk\p(t)\Xk ) , (77) 

k 

where {\Xk)}k is any basis of Hp. 

For this reason, we will represent each site with a two- 
dimensional Hilbert space H& (k = 1,2, 3,4) with basis 
{| rife = 0), | rik = 1)} indicating the presence (n^ = 1) or 
absence (n^ = 0) of a random walker. Since |V^) denotes 
a random walker at site fc, we can write the site basis as 

IVh) = |m = 1) 0 |n 2 = 0) 0 |n 3 = 0) 0 |n 4 = 0) , (78) 

1^2) = |m = 0) 0 |n 2 = 1) 0 |n 3 = 0) 0 |n 4 = 0) , (79) 

1 ^ 3 ) = |ni = 0) 0 |n 2 = 0) 0 |n 3 = 1) ( 8 ) |n 4 = 0) , (80) 


p (3) (£ + dt) = (l - k 2 i dt) |'0i)('0i| . 
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IV^) = \ni = 0) 0 \n 2 = 0) 0 |n 3 = 0) 0 |n 4 = 1) . (81) 


This method of defining sites can be viewed as defining a 
four-mode state in quantum optics with each mode con¬ 
taining at most one photon (or equivalently a four-mode 
fermionic state). Note that this expresses M in the form 


of (76) with 


Mr = Mi 0 M 3 , Mp = M 2 0 M 4 . (82) 


1 1 

= (ri2 = m, 714 = m'\p(t)\n 2 = m, 774 = m!) 

m =0 m '=0 

(84) 

= (n 2 = 0, n 4 = 0| p{t) | n 2 = 0, n 4 = 0) 

+ (n 2 = 0 , n 4 = 1| pit) \n 2 = 0 , n 4 = 1) 

+ (n 2 = 1 , n 4 = 0| p(t) |n 2 = 1 , n 4 = 0) . (85) 


It will be convenient to introduce a short-hand notation 
in which the tensor product is omitted. In general, we 
will write 


\nj,n k ) = \nj) ® \n k ) . (83) 

The partial trace of ( |49| ) over sites 2 and 4 is then given 

by 

PR(t) 


We have noted the expectation of p(t) with respect to the 
state | ri 2 = 1,714 = 1) is identically zero since there can 
be at most one particle (or walker) in the system. The 
calculation of (85) is a little bit involved so we leave the 
details to Appendix [B] The result is however simple to 
state. Noting that /?r now refers only to 72 3 )} nijn3 , 
we will simply denote a state with n% = 0 and 77,3 = 1 as 
|0,1) and write 


Pr(£) — £qw Pr(£) = [&21 Pn(t) + /C 43 P 33 (t)~\ |0,0)(0,0| — A : 43 P 3 3 (t) |0,1)(0,1| — k 2 1 pn(t) |1,0)(1, 0| 
- I (k 21 + fc 43 ) p 13 (t) |1,0)(0,1| - I (k 2 1 + ^ 43 ) P3i(0 |0,1)(1,0| . 

As before, this equation is perhaps easier to read from its matrix representation, given by 
(k 2 i pn + k 2 3 p 33 0 0 


PR = 


V 


0 

0 — k 2 spss — §(^21 + ^23)^31 | = 

0 — | (&21 + ^23) P13 — k 2 \ pn 


( 86 ) 


(0,0|p R |0,0) (0,0 |pr| 0, 1) (0,0|p R |l,0) 

(0,1|Pr|0,0) (0,1|/9r|0,1) (0,1|pr|1,0) 

(1> 0|pr|0, 0) (1, 0|pr|0, 1) (1, 0|pr|1, 0) 


(87) 


This gives the correct rates for dephasing and popula¬ 
tion transfer for the radical pair except now we have one 
additional state, |0,0), whose effect is to drain the pop¬ 
ulations out of |0,1) and |1,0). It would be tempting to 
compare the dynamics given by ( |87| ) to the graph of Fig. [4] 


and define |S) = |0,1), |T) = |1,0), and |P) = |0,0). 
However, the partial trace does not correspond to this 
identification, because as we have just shown above, the 
evolution defined by Fig. [4] is given by (59) [or equiva¬ 
lently by (60)—(68)] which has the matrix representation 


P = 


^ &21 Pll + ^23 P33 ~ \ &23 P23 ~~ \ &21 P21 

^23 P33 — \ {k 2 \ + k 2 %) P 31 

— \ (k 2 \ + k 2 z) P 13 — k 2 \ pn 


— \ k 2 3 P32 

V —\k 21 p 12 


<P|p|P) <P|p|T> <P|p|S) 
<T|p|P) (T\p\T) <T|p|S) 
\<S|p|P) (S\p\T) (S|p|S) 


( 88 ) 


r 


Note that we have used (43) and (57) and reordered states whereas (87) does not. The two matrices do not 

the matrix elements for ease of comparison with (87). even refer to the same basis—we are correct to equate 

The difference between (87) and (88) is clear. Equation 11,0) to |S), and |0,1) to |T) in ( |87| ), but we would be 

(88) has coherences between the product and radical-pair mistaken to identify |0, 0) with |P). The reason is because 
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Fig. 5: Depiction of the different approaches to the radical-pair reaction operator based on what states are included in the model. Each 
box is colour coded in accordance to the colour coding of Figs. [3] and [4] Each site will have either a filled dot or an unfilled dot, representing 
respectively the presence or absence of the random walker, (a) The conventional/Haberkorn model with the reaction operator referring to 
only the singlet and triplet states of the radical pair. This does not preserve the norm of p{t) and can be derived from (b) by discarding 
|P S ) and |Pt)- (b) Including the reaction products in the reaction operator so that p{t) is normalised at all times, (c) The partial-trace 
method which expresses our ignorance of the products. This gives rise to a state in which the radical pair is neither the singlet nor triplet. 
We represent the ignorance of products by question marks at the product sites. 


|0, 0) conveys no other information except the absence of 
the random walker from sites one and three. It does not 
say where the walker is. Thus |0,0) should be regarded 
as a radical-pair state because it gives us only informa¬ 
tion about the radical pair, namely that it is in neither 
the singlet nor triplet state, consistent with the fact that 
we have traced over the products in (88). In contrast, 
|P) says exactly which site the random walker is at. It 
is represented by a node on the graph in Fig. [4] whereas 
|0,0) is not. For notational consistency we will write 


|N) = |0,0) , (89) 

where N may stand for “neither”, “none”, or “null”. 
This distincition between |0,0) and |P) is an impor¬ 
tant and interesting one because it suggests that |N) 
is another radical-pair state that we should consider 
and therefore extend the minimal basis from {|S), |T)} 
to {|S), |T), |N)}. Previous treatments on the radical- 
pair reaction operator have been to use either a trace- 
decreasing p{t) without products, or a trace-preserving 
p{t) with products. The partial trace has the advan¬ 
tage that it is both trace-preserving and excludes the 


products. It achieves this by regarding |N) as just an¬ 
other radical-pair state which has not been considered 
(or taken seriously) before. We summarise the previous 
approaches to radical-pair kinetics alongside the partial- 
trace method schematically in Fig. [5j 

Finally, we make a couple of observations of the model 
defined by (86) and (87): 


1. We first note some similarities and differences be¬ 
tween the partial-trace approach and the conven¬ 
tional Haberkorn model. The two models are es¬ 
sentially equivalent in that the former description 
uses a 3 x 3 matrix with unit trace, whereas the 
latter uses a 2 x 2 matrix plus one scalar equation: 


pss{t) + Ptt(£) = —ks f dt'pss{t') — &t [ dt'pTT{t') • 
Jo Jo 

(90) 


We have used (51) to express the product popula¬ 


tions on the right-hand side in terms of pss{t) and 
Ptt(^)- However, the right-hand side of ( |9Q| ) explic¬ 
itly refers to product populations whereas the trace 
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of does not. This is because |N) is a radical-pair 
state so that Tt[pr] is interpreted as the probabil¬ 
ity of finding the radical pair to be either a singlet, 


triplet, or neither. Equation (90) provides more 


detail by saying the probability of finding the sys¬ 
tem in either the singlet, triplet, singlet product, 
or triplet product states must be one. Naturally, 
the probability of not being in the singlet or triplet 
state must be equal to the probability of being in 
the product states. 

2. We mentioned that the random walk of Fig. [4] is 
not the same as the partial trace because (|87|) and 


(88) are not equivalent. However, if we note that 


coherences between the radical pair and products 


in (88) are decoupled from the rest of the matrix, 


then the random walk of Fig. [4] can be used to sim¬ 
ulate the partial trace by first defining |^i) = |S), 
|^ 2 ) = |N), and 1 ^ 3 ) = |T), and then setting 
(A I p(t)\fo) = (^3 |p W 1^2) = 0 in the end. This 
defines the partial trace over products in terms of 
a random walk. 


V. EXTENSION TO MORE GENERAL CASES 


H 31 

A- 


IV’i) • < 


1 


=> • IV>3> 


721 


£*31 


743 


1 ^ 2 ) 


IV>4> 


Fig. 6 : A quantum walk involving all three processes described 
in Sec. |III| Coherent evolution is represented by a green two-way 
arrow and parameterised by £31 (recall that the actual rate of os¬ 
cillation is given by 0 : 31 ). This models the effect of the magnetic 
interactions such as the Zeeman or hyperfine. Any extra dephasing 
(i.e. not caused by recombination) can be effectively described by 
the dephasing map. If applied to the experiment of Ref. m this 
would model the ^-anisotropy for one of the radicals. 


by the recombination processes alone. Although the ex¬ 
periment minimises as many decoherent processes as pos¬ 
sible not all such processes can be made negligible. We 
therefore include an additional source of dephasing in our 
model. This gives the resultant graph in Fig. [6] whose 
time evolution can be described by the map 


Not all processes in a chemical reaction can be rep¬ 
resented by amplitude damping. In this section we con¬ 
sider reactions where coherent oscillations and additional 
dephasing can occur. This is the case with the radical- 
pair reaction where magnetic interactions give rise to co¬ 
herent oscillations between the singlet and triplet states. 


It should be clear that (46)—(48) can be generalised to 


any concatenation of maps describing amplitude damp¬ 
ing, dephasing, or unitary evolution. More precisely, a 
graph defined by 

= {Ugh ' ' • U c d U a b) ( V xy • • • V rs Vpq ) 

X (M mn Mu Mu) , (91) 

can be written as /C(t) = exp(^t) with 

Q — (j^gh + ’ ’ ’ + Ttcd + T^ab) + {S X y + ' ' ' + $rs + Spq) 

+ {£mn + ' ' ' + AcZ + Aj) • (92) 


We have omitted the time argument in (91) for clarity 


and used (32), (33), (38), and (39) to obtain (92) 


A. Application to experiment 


The experiment of S ec. Ill B 2 is an example where a 


model in the form of (91) and (92) would be directly 
applicable. Recall from Sec. |IIB 2 that the measured de¬ 


^Qw(dt) = V3i{dt)Usi(dt) A44s(dt) M.2i{dt) . (93) 

Applying ( [91] ) and ( [92] ) gives the reaction operator: 

Gqw p(t) = 


+ &21 
+ &43 
+ 431 


Q21 p(t) Qh - 1 Qi p(t) - 1 p(t) Q\ 
Q43 p(t) q \ 3 ~\Qs pit ) - \ pit) Qs 
Qi pit ) Qi- | Qi pit ) - 1 pit) Qi 

(94) 


The last line in (94) models any additional dephasing oc¬ 


curring on top of the recombination without contribut¬ 
ing to population transfer (recall Sec. |IIIB| ). Maeda and 
coworkers have suggested such a process [31]. An order- 
of-magnitude estimate suggests the additional dephasing 
observed in their experiment can be accounted for by the 
anisotropy associated with the electron gyromagnetic ra¬ 
tio for one of the radicals (also referred to as ^-anisotropy 
in electron paramagnetic resonance) m- Their order- 
of-magnitude estimate for the rate of dephasing due to 


^-anisotropy could in principle be used in (94) for q%i. 


We will study of the effect of a nonzero q%i in a sequel 
paper where a toy model is used to simulate radical-pair 
reactions in plant cryptochromes m- However, there 
we allow (731 to be a free parameter so the dependece of 
the reaction kinetics on additional dephasing can be seen 


phasing rate does not correspond to the dephasing caused 
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We noted earlier in Sec. IIII Al that the order of sub¬ 
scripts in the amplitude damping map is important as 
this determines the direction of the process. This is also 
the case for the dephasing map and is obviously true 


from (32) as its generator is determined solely by the 


projector corresponding to the first index (counting from 
right to left). If we include products into our reaction 


operator then the form of (32) will not only affect the 


singlet-triplet coherence of the radical pair, but also the 
coherences between the radical pair and products. The 
ordering of the indices then determines which coherences 
between the radical-pair and product are affected. Take 
the four-state graph in Fig. [ 6 ] for example, the contribu¬ 
tion to the reaction operator given by V 31 (dt) is 


<?31(#31 )P 

( 0 


“2 #31 P 21 
~ \ #3l/>31 

V4 931 P41 


"2 #31 Pl2 
0 
0 
0 


- \ # 3 lPl 3 — \ #31 P 14 


(95) 


whereas the contribution from V 13 (dt) is 

*^13(#13) P 

( 


0 

0 

- \ #13 P31 

0 


0 

0 

"2 #13 P32 

0 


" 2 #13 Pl3 
- \ #13 P23 
0 

‘ 2 #13 P43 


0 

0 

- 2 #13 P34 

0 


(96) 


Note that we have explicitly written out the dependence 
on the dephasing rate in the generators to be clear. This 
therefore raises the question as to whether a dephasing 
map with a particular ordering of indices should be pre¬ 
ferred over another, e.g. V 31 (dt) as opposed to V 13 (dt) in 
(94). It turns out that this asymmetry is usually not a 


problem for describing radical-pair reactions because the 
radical pair is usually created in either the singlet or the 
triplet state. More generally, the asymmetry in the de¬ 
phasing map is irrelevant for any initially mixed state of 
the form 


N 


p( 0 ) = ^2pk iV’fcXV’fci 


(97) 


k =1 


This is the case with the experiment of Sec. IIB 2 where 
the radical pair begins in the singlet state, which corre¬ 


sponds to (97) with N = 4, pi = 1 , and P 2 = p 3 — 
p 4 = 0. This of course assumes that there are no other 
processes present which create coherences between the 
radical pair and products. We comment further on the 
asymmetry of the dephasing map and its use in deriving 
reaction operators below [56]. 


P31 P'13 

AV 

• • 

l^i) 1^3) 


Fig. 7: We represent two applications of the dephasing map by 
using multiple triangles with a single line. This gives an overall 
map that is symmetric in the decoherence rates when the rates for 
each direction are set equal. This is also the Kominis model with 
the dephasing rates replaced by the recombination rates. 


B. Variant models 


1. Symmetric dephasing 


If one insists that a symmetric dephasing model be 
used then this can be accomplished by two applications 
of the dephasing map. We represent such an operation 
by Fig. Ill Since we associate an upwards-pointing tri¬ 
angle with the map Vsi(qsi m ,dt) it makes sense to use 
a downwards-pointing triangle for Vis(qis;dt). We thus 
define the “direction” of dephasing by the direction in 
which the triangle points. A symmetric dephasing map 
may then be defined by 


W 31 (q;dt) = V 3 i(<? 3 i := q;dt)Vi 3 (qi3 ■= q\dt) . (98) 


The generator corresponding to this is simply given by 
the sum of ( |95| ) and (J96| with q^i = #13 = #: 

A 31 (#) p = [<Si 3 (#) + S 3 i(#)]p 


/ 

0 

-\qp12 

-qpi 3 

-\qpu\ 


-\qp21 

0 

-\qp23 

0 


-qp 3 i 

-\qp32 

0 

-\qp 2 A 

V 

-\qpu 

0 

-\qpv> 

0 / 


(99) 


We can also provide a model of dephasing in which only 
the radical pair is referred to by using the partial trace 
defined earlier in Sec. |IVB 2 [ Since the partial trace of a 
sum is the sum of partial traces we can simply take the 


partial trace of (99) and add it to the generator for other 


processes. We can use exactly the same procedure as 


before to calculate the partial trace of (99) but it should 
be intuitive that the result is given by 


Pr = 


too 0 

0 0 -qp3i 

\0 ~qpi3 0 


( 100 ) 


where ( 100 ) is expressed in the {|ni, ns)} ni ,n 3 basis and 


the ordering of the matrix elements is the same as (87). 
Adding ( | 100 ) to ( |87| ) thus gives a model for the recombi¬ 
nation process with added dephasing in the radical-pair 
basis {|S),|T),|N)}. 
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2. Relation to Kominis’s reaction operator 

This asymmetry in the dephasing map can be used 
to obtain Kominis’s model operationally. Recall from 
Sec. IIB1 that his reaction operator is given by @ 
and is trace-preserving, with the radical-pair population 
evolved using a separate equation, given by (§. This 
model can be seen as two applications of the dephasing 
map in opposite directions: 

p(t + dt ) = V3i(<?3i := ks; dt) Vi 3 (gi3 := k T ; dt) p(t) . 

( 101 ) 


Using (43) this is equivalent to 


ks) <Si3((Zi3 &t) p(t) 


= k s 
+ Ajt 


Qs p(t) Qs-^Qs p(t ) - ^ pit) Qs 


Qt Pit) Qt — 2 Qt Pit) - - pit) Qt 


= C-k pit) • 


( 102 ) 


We can now understand why Kominis had to introduce 
an ad hoc method for describing population loss in the 


radical-pair reaction, namely that (102) only describes 


the loss of coherences. Of course we have not analysed 
the physics of Kominis’s model. All we have done is 
to point out that whatever physics and assumptions go 
into Kominis’s model, they amount to dephasing for the 
radical pair. This should be compared with the quantum- 
walk model used here in which the radical-pair dephasing 
is seen explicitly as a consequence of the decay out of |Vq) 
and ^3) to products. If we want to describe population 
loss in the minimal basis {| S), |T)} then it makes more 
sense to use the conventional model given by Cup (or 
equivalently Cqyjp). We can of course also tack on Cup 
to (102) to arrive at a model where there is population 


loss and additional dephasing in the minimal basis: 


dp - (r 


£h) pit) 


(103) 


The graph corresponding to Kominis’s model is also given 
by Fig. 13 when the dephasing rates are set equal to the 
recombination rates. 


VI. SUMMARY AND DISCUSSION 

We have shown that quantum walks provide the same 
level of treatment for coherent chemical kinetics as rate- 
equation models do for classical chemical kinetics by ap¬ 
plying it to the radical-pair reaction of the avian com¬ 
pass. However the quantum walk considered here is not 


the same as those in the quantum-walk literature where 
an additional system degree of freedom is used as a coin 
Email. The simplicity of our approach lies in the de¬ 
composition a multisite reaction into two-site processes. 
If multiple processes occur between two sites then we 
can also consider these processes separately. This gives 
a systematic method of deriving reaction operators. The 
breakdown into two-site processes also makes our ap¬ 
proach less prone to modelling errors since two-site pro¬ 
cesses are much simpler to study—they are in fact well 
known in quantum information theory. As we saw in 
Sec. IV| this allowed us to obtain a recombination de¬ 
phasing rate for the standard radical-pair reaction which 
is consistent with experiments. This follows a similar 
line of thought as the Jones-Hore work except they at¬ 


tempted a derivation “in one go”. Section IV can also be 
seen as a derivation of the Haberkorn reaction operator 
using a theory of coherent chemical kinetics as argued. 
We have also shown how the partial trace can be used 
to obtain a model where the products are ignored but 
still preserves the normalisation of the radical-pair state. 
This gives rise to a third state that has not been men¬ 
tioned before in the radical-pair literature and can be 
seen as a half-way approach to the existing models (see 
Fig# We have also considered more general reactions 
in Sec. [V] where a model corresponding to the experiment 
of Ref. [31] was introduced. We also considered variant 
forms of dephasing and discussed its relevance to the Ko¬ 
minis model. 

Due to space we have not provided the details of how 


additional dephasing like the one in the model of (94) will 


affect a coherent reaction. The effect of such decoherent 
processes will ultimately depend on the actual values of 
the intermediate transition rates. In a sequel paper we 
further illustrate the use of quantum walks by simulat¬ 
ing a toy model for a radical-pair reaction in plant cryp¬ 
tochromes where order-of-magnitude estimates for the in¬ 
termediate transition rates are available m- However, 
instead of using a particular value of dephasing, we will 
vary the dephasing parameter over its full range so as 
to study the deviation of a maximally coherent reaction 
from one that is fully classical. For this purpose we use 
a well-known quantity called the hitting time from quan¬ 
tum walks as a measure of the quantumness of the reac¬ 
tion. This measures how quickly the reaction occurs as 
a function of the coherences in the reaction. Our sequel 
paper thus shows how the “full” machinery of quantum 
walks can be applied, not just the idea of composing maps 
presented in this paper. Such an application of coherent 
chemical kinetics can then be used to argue if a quantum 
description for the reaction being studied is necessary, or 
if simple classical rate equations will do. 

Aside from the Jones-Hore work, the open-systems ap¬ 
proach due to Tiersch and coworkers [29 also deserves 
some comment. There is a sense in which our paper is 
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similar to their’s and a sense in which it is the complete 
opposite. Our work is similar to Ref. [29] in that Tiersch 
and colleagues also use maps to describe different sorts of 
evolution for the system (i.e. the radical pair). Once the 
physics is seen clearly the system evolution is then de¬ 
scribed by a master equation, derived by converting the 
maps into differential form. However, whereas our goal 
is to establish the systematic use of maps for deriving 
reaction operators because of their operational nature, 
Tiersch and others use maps to suggest what underly¬ 
ing physics might give rise to known reaction operators. 
More precisely, Tiersch and colleagues treat the radical 
pair as an open system and ask what sort of system- 
environment interactions can give rise to sensible reac¬ 
tion operators. Their maps are thus different to our’s 
in that their’s act on the joint space of the system plus 
environment. In contrast, our maps act directly on the 
system alone. One can therefore say that our work and 
Ref. [29] are similar in methodology due to the simplicity 
of maps, but quite different in spirit. 
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Appendices 

Appendix A: Amplitude-damping Kraus map 
1. Operator-sum form 

For ease of reference we repeat the amplitude-damping 
map here: 

p(t + At) = M 2 i p{t ) 

= M$iAt)pit)[M$iAt )] f 

+ M^iAt)pit)[M^iAt)}\ (Al) 

where the two Kraus operators are 

V21 \ At ) = Vl 2 i(M) IV^XV’iI , (A2) 

M^fiAt) = |- 02 ><- 02 | + v / l-72i(At) \ipi)(ipi\ ■ (A3) 


The forms of M^iAt) and M^iAt) can be understood 
by first noting that during At we either observe a tran¬ 
sition or we do not (assuming of course that the system 
can be probed to give us this information). Considering 
first the case of observing a transition, the state is then 

p {1) it + At) = (At) pit)[M^ (At)] 1 /pi (At) . 

(A4) 

The probability of observing the transition is simply 
given by 

Pi (At) = 72i(Af)pii(f) , (A5) 

where pjj(t ) = (^j\p(t)\^j) is the probability of finding 
the system to be in \^j) (j = 1,2) at time t. It should be 
clear directly from the form of (At) that it effects 
the required conditional change as it is proportional to 
(VhI? whose effect is to take a system in state |^i) to 
the state IV^)- The inclusion of ^721 (At) then allows 
the correct probability of observing a transition to be 
found according to the formalism. Note the expression 
for pi (At) is exactly what one would expect from clas¬ 
sical probability theory. Considering now the case when 
we do not observe a transition, the state is then given by 

p {2) it + At) = (At) pit)[M$ ( A U / p2(Af) , 

(A6) 

where 

P 2 ( At) = P 22 it) + [1 - 721 (Af)] pn it) , (A7) 

is the probability of not observing a transition. These 
results can be understood by noting that two possible 
scenarios contribute to a no-transition observation: Ei¬ 
ther the system is already in state ^2} at time t, which 
occurs with probability P 22 (f)•> or it is h 1 state \^\) but 
has not yet jumped to 1^2), which occurs with proba¬ 
bility [1 — 72i(At)]pn(£). The probability of not seeing 
a transition is therefore the sum of the probabilities for 
each of these scenarios. That (At) describes a com¬ 
bination of these two scenarios can also be seen from 
its form, which we can understand by a simple analogy 
to M^(At). Consider first the case where the system 
is in state |^i) and remains in |Vq)- Instead of taking 
|^i) to |^2) as in M^(At), we now take |^i) to itself. 
This means that we simply replace the transition opera¬ 
tor 1^2)(Vh| in M^(At) by the projector |Vh)(Vh|- We 
would also have to replace the 721 (At) under the square 
root in M^\At) by 1 — 721 (At) since now we are con¬ 
cerned with the case where the system stays in Do¬ 
ing so gives us M^ 2 )(At) = yjl — 721 (At) |Vh)(Vh|, but 
we know this is not the complete description yet as we 
have not considered the contribution due to the system 
being in 1^2) and staying there. If the system is already 
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in |^ 2 ) 5 the process should take ^ 2 ) to itself because no 
other processes are present which can take the system 
out of l^)- The probability that the system remains in 
^ 2 ) given that it was in ^ 2 ) is thus simply 1. Therefore 
we simply add the projector ^ 2 ) (V^l (with coefficient 1 ) 
to a/T — 721 (At) |' 0 i)(V ; i| to arrive at the resultant form 
of (At). It is trivial to show that M%i (At) produces 
the correct state by letting p(t) be |Vh)(Vh| and 1^2)(V^l 
in turn. 

For the mathematically inclined reader, we note that 


(A3) can be derived directly by using (A 2 ) and the con¬ 
straint 

[Mg (At)]'Mg (At) + [Mg (At)]'Mg (At) = 1 . 

(A 8 ) 


By resolving the identity on the right-hand side in the 
site basis, ( |A 8 | ) gives 

[Mg (At)]'Mg (At) = [l -72i(Ai)] |V’i)(V’i| 

+ |V>2> (lp2 1 • (A9) 


Note that Mg (At) is simply the operator square root of 
this equation. Since ( |A9| ) is diagonal, we arrive at ( |A3| ) 
on taking the sqaure root of the coefficients of |Vq)(Vq| 
and |^ 2 )(* 02 1- This derivation of (A3) is simple but does 
contain the insight provided above. 

Now that we have the necessary Kraus operators the 
evolution of the system follows directly by forming the 
sum ( |A 1 | ). Notice from the above that (At) and 

/ON 

(At) are essentially time-evolution operators but 
conditioned on our knowledge of whether the system un¬ 


derwent a jump or not. Equations (A4) and (A 6 ) are 


in fact quantum analogues of the classical Bayes rule 
[53] , Conditioning requires that we monitor the system 
for the entire duration of At. What the operator-sum 
representation of p(t + At) describes is how the state 
should evolve without us having to monitor the system 
continuously, or in the language of probability theory, 
it describes the unconditioned state. This can be un¬ 
derstood as follows: If one does not monitor the system 
then all we can say is that with probability pi (At) the 
system will be in the state p^\t + At), and with proba¬ 
bility p 2 (At) the system will be in the state p( 2 \t + At). 
From probability theory, we would say that the state in 
the absence of such monitoring at time t + At (i.e. the 
unconditioned state) is therefore a weighted sum of the 
conditioned states p^(t + At) and p^ 2 \t + At), 

p(t + At) = pi (At) p^\t + At) + p 2 (At) p( 2 \t + At) . 

(A 10 ) 


Substituting in (A4) and (A 6 ) we obtain exactly the re¬ 


sult of Kraus. In practice one often has an ensemble 
of particles and all we know is the fraction of particles 
that underwent a state transition during At. In this case 


pi (At) is simply the fraction of particles that jumped 
and p 2 (At) the fraction that did not. Note also the dif¬ 
ference between pi (At) and 721 (At)—The former is the 
probability of observing a transition from t to t + At 
without assuming that we know which state the system 
is in at time t, whereas the latter does, 721 (At) is the 
probability of jumping conditioned on the system being 
in state 17 / 7 ) at time t. 


2. Differential form 


We have now described a simple one-way population 
transfer completely as a probabilistic process. Instead 
of expressing the system evolution as a sum over condi¬ 
tioned states we can express it in the form of a differential 
equation. Such an equation can be derived by consider¬ 
ing the evolution of p(t) over an infinitesimal interval dt. 
In this case it is more appropriate to refer to the rate 
at which the system jumps from |^i) to |^ 2 ) over some 
interval At rather than the probability 721 (At). If we de¬ 
note the fraction of particles that jump from |^i) to ^ 2 ) 
per second by & 2 i> the probability 721 (At) is then related 
to &21 by 721 (At) = Jfe 2 iAt. This means that the proba¬ 
bility of a jump in an infinitesimally small time interval 
is also an infinitesimal. The evolution of the system state 
is now given by 

p(t + dt) = Mg (dt) p{t) [Mg (dt)]' 

+ Mg(dt)p(t)[Mg(dt)]' , (All) 

Using the binomial expansion we have 

\A - 721 (dt) = l--^dt. (A 12 ) 

We can then write the Kraus operators for infinitesimal 
evolution using Q 21 = IV^)(Vh I and Q 1 = IVhXVhl as 


M^(dt) = \Jk 2 i dt Q21 , (A13) 

Mg (dt) = i -^yQidt. (A14) 


Substituting (A13) and (A14) into (All) and neglecting 
terms on the order of dt 2 we arrive at 


dp 

dt 


— &21 


Q21 p(t) q ' 21 


\ Qi p(t) ~ \ p(t) Qi ■ 

(A15) 


This is the master equation corresponding to the am¬ 
plitude damping map and can be put in the Lindblad 
form if one wishes by using the property Q 1 = 1 Q 2 i- 

We have derived this equation rather simply by apply¬ 
ing the Kraus formalism, hence it can be regarded as 
merely a restatement of the operator-sum representation 
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of p(t) in differential form. It is only a matter of pref¬ 
erence whether one wants to use a map or a differential 
equation to simulate the system dynamics but the Kraus 
formalism provides a simple way to understand the essen¬ 
tial physics of the process by using only basic probability 
ideas. 


Appendix B: Partial trace over chemical products 


Here we show how to obtain an equation of motion for 
the radical pair where the products are ignored by using 
the partial trace operation on (49). We have already 
argued in Sec. |IVB~2| that this is given by 


hi{t) = {n 2 = 0, n 4 = 0| p(t) | n 2 = 0, n 4 = 0) 

+ (n 2 = 0,n 4 = 1| p(t) | n 2 = 0,n 4 = 1) 

+ (n 2 = 1, n 4 = 0| p(t) | n 2 = 1, n 4 = 0) . (Bl) 


where 


p(t) = k 2 1 

+ &43 


Q21 p{t) Q f 21 * I Qi p{t) - ^ pit) Qx 


<?43 Pit) Q \ 3 - - Qz pit) - - pit) Q 3 


(B2) 


and we have written the site basis using 177 - 1 , 77 - 2 , 77 - 3 , 77 - 4 ) 
as [recall (J78|-(81) and ( |83| )] 

IVq) = 11,0,0,0), 1 ^ 3 ) = 10,0,1,0), (B3) 

1 ^ 2 ) = 10,1,0,0), 1 ^ 4 ) = 10,0,0,1). (B4) 

The sum ( |B1| ) is most easily calculated by first rewrit¬ 
ing (B2) as 

/ X 4 

P = k 21 i p U 1^2) (^2 | - -Z Plm |V’l)(V’m| 

\ m= 1 

- \ ^2 P ml I i>m)(i’l\^ + ^43 ^33 ^ 4 ) (^41 

1 4 1 4 \ 

- P3m IV’3)(V’m| ~ ^^2 Pm3 | V’mXV’sl • 


m= 1 


m= 1 


(B5) 


— $ m2 (1 ^777,4) 177-1 — — ^m3) • (B8) 

We will also use the summation convention where a re¬ 
peated index is summed over. Using (B5) the first term 
in 


can thus be calculated as follows. 

(n 2 = 0, n 4 = 0| p(t) | n 2 = 0, 77-4 = 0) 

1 


~Y (pim («2 = 0,n 4 = 0 \tpi) (tp m \n 2 = 0 ,n 4 = 0) 
+ Pmi (n 2 = 0, n 4 = 0 1 ip m ) {ipi\n 2 = 1, n 4 = 0) J 
- ^ (p3m («2 = 0,n 4 = 0|^ 3 ) ii>m\n 2 = 0 ,n 4 = 0) 
+ Pm3{n 2 = 0,n 4 = 0|t/> m )(i/> 3 |n 2 = 0,n 4 = 0)) , 


(B9) 


where we have noted on using (B6) that 
(n 2 = 0, n 4 = 0 |V^ 2 ) = (ri 2 = 0, n 4 = 0 |^ 4 ) = 0 . (B10) 
while the first term in (|B9|) is 


Plm (n 2 = 0 ,n 4 = 0 \ipi) (ip m \n 2 = 0 ,n 4 = 0 ) 
= Plm (1 ) (1 5 m A ) 

x |m = 1 ,n 3 = 0 )(ni = 5 m i,n 3 = < 5 m3 | 


— Plm (1 5 rn 4 5 m2 llmT^rnA) 

x |m = l,n 3 = 0)(ni = 5 m i,n 3 = 5 m3 | 


= |m = l,n 3 = 0) yPlm {n% = 5 m i, n 3 = 6 m3 1 
5mA Plm (ill ^ml 5 H3 5m3 1 

5 m2 plm (ill = 5rn 1 ,II 3 = 5 m 3 1 J 


= Im = 1, n 3 = 0) (pii (rn = 1, n 3 = 0| 



+ P 13 (hi = 0 , n 3 


1| 




= P11 |ni = 1,113 = 0)(ni = l,n 3 = 0| 

- P 13 Im = 1 , 113 = 0)(m = 0. n 3 = 1 . (Bll) 


This also gives us the second term in (B9) since it is just 
the Hermitian conjugate of ( Bll| ), 


The following identities will therefore prove useful 


(112 = 

0 , 77-4 — 0 I'lprn) 



= (1 

— ^m2) (1 — ^m 4 ) 

1 77-1 — ^777,1 ? 77-3 — ^ 777 , 3 ) 5 

(B 6 ) 

(n 2 = 

0 , 77-4 — 1 | 'Ipm) 



= (1 

$rri2 ) ^ 777,4 |t7-i — 

' ^ml,77-3 — ^ 7773 ) 5 

(B7) 

(n 2 = 

1 , 77-4 = 0 1 2pm) 




Pm 1 (n 2 = 0, n 4 = 0 |?/> m ) (f/>i |n 2 = 0, n 4 = 0) 

= (pirn (H2 = 0,n 4 = 0 |^i)(^m|n 2 = 0,n 4 = 0)) 

= Pn |ni = 1,H3 = 0)(ni = 1,n 3 = 0| 

— P 3 i Im = 0, n 3 = l)(iii = 1,n 3 = 0| . (B12) 

The third term in ( |B9| ) can be calculated in exactly the 
same manner with p^ m replacing pi m , and I'lps) replacing 
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I'lpi). By inspection of the above working we see that this 
amounts to making the following replacements in (Bll) 
(bras remain unchanged): 


Pi 1 —> P31 5 Pi3 —> P33 5 (B13) 

|ni = l,n 3 = 0) —» \rii = 0,n 3 = 1) . (B14) 

We thus obtain 

P 3 m (n 2 = 0,n 4 = 0 1 ^ 3 ) (V’rn |n 2 = 0,rc 4 = 0) 

= P 31 = 0,n 3 = l)(ni = l,n 3 = 0| 

- P 33 K = 0,n 3 = l)(ni = 0,n 3 = 1| , (B15) 


and taking the Hermitian conjugate, 


Pm 3 (n 2 = 0, n 4 = 0 \'lprn) (^ 3^2 = 0, n 4 = 0) 

= (p 3 m (n 2 = 0, n 4 = 0 |V^ 3 ) (^rn|n 2 = 0, n 4 = 0)) 


= P 13 \ni = 1 ,n 3 = 0)(ni = 0,n 3 = 1| 

- P 33 |ni = 0,n 3 = l)(ni = 0,n 3 = 1| . (B16) 


Substituting ( |B11 ), ( B12[ ), ( |B15[ ), and (B16) into ( |B9| ) 
and collecting like terms we arrive at 


(n 2 = 0, n 4 = 0| p(£) |n 2 = 0, n 4 = 0) 

= ~k 2 ipii |ni = l,n 3 = 0)(ni = l,n 3 = 0| 

- ^ (fc 2 i + & 43 )pi 3 K = l,n 3 = 0)(ni = 0,n 3 = 1| 

- ^ (k 2 1 + A: 43 )p3i = 0,n 3 = l)(ni = l,n 3 = 0| 

- ^43 P 33 Wi = 0, n 3 = l)(ni = 0, n 3 = 1| . (B17) 


The remaining two terms in ( p5[ ) are much easier to 
calculate with the help of (B7) and (B8). The second 
term in (|B1|) is given by 


(n 2 = 0, n 4 = 11 p(t) \n 2 = 0, n 4 = 1) 


— k^z P 33 \ni — 0,n 3 — 0)(ni — 0,n 3 — 0| , (B18) 


where we have noted from (B7) that 


(n 2 = 0, n 4 = 1 |^m) = |m = 0, n 3 = 0) for m = 4 , 

= 0 form = 1,2, 3. (B19) 


and 


(n 2 = 0, n 4 = 11^4) = |ni = 0, n 3 = 0) . (B20) 


Similarly, the third term in (Bl) is 


(n 2 = 1, n 4 = 0| p(t) | n 2 = 1, n 4 = 0) 

= AJ 21 P 11 = 0,n 3 = 0)(ni = 0,n 3 = 0| , (B21) 


where we have noted from (B8) that 


(n 2 = 1, n 4 = 0 IV’m) = |ni = 0, n 3 = 0) for m = 2 , 

= 0 for m = 1,3,4. (B22) 


Collecting (B17), (B18), and (B21), we finally have the 
state pr which has the products traced out. The final 
form of (Bl) is then 


Pr = (k 21 pn + ^43 P 33 ) |0,0) (0, 0| - & 43 P 33 |0,1)(0,1| 

- k 2 1 Pll |1,0)(1,0| - i (k 2 1 + fc 43 ) pis |1,0)(0,1| 
-1 (fc2i+fc 43 )p3i|0,l)<l,0| . (B23) 

where we have used the fact that pr must be spanned by 
{|^i 5 ^3)}ni,n 3 to omit writing out ni and n 3 explicitly 
in ( |B23[ ). 
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